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Abstract 



The electron-electron relaxation in impure two-dimensional superconduc- 
tors is studied. All channels of the electron-electron interaction classified in 
the Nambu representation are taken into account. It is shown that the re- 
combination relaxation rate originates from quasipartical processes associated 
with fluctuations of the electron density and the phase of the order parameter. 
At low temperatures the recombination relaxation rate has a double exponen- 
tial temperature dependence. The scattering relaxation rate at low temper- 
atures has a power law temperature dependence due to contributions from 
gapless collective excitations, the phase modes. Two-layer superconductor- 
normal metal system is also considered. It is shown that the recombination 
relaxation rate in the superconducting layer has a single exponential factor 
at low temperatures in comparison with a one layer superconducting system. 
This increase in the recombination relaxation rate originates from the inter- 
layer Coulomb interaction and may be used in constructing of superconducting 
radiation detectors. 

PACS: 73.50.Bk, 73.50.Dn. 
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I. INTRODUCTION 



The electron-electron energy-relaxation time determines a number of parameters of 
nonequilibrium superconductors such as the relaxation times for the amplitude and the 
phase of the order parameter It is also important for nonequilibrium superconducting 
radiation detectors based on the resistive and inductive responses. The electron-electron 
relaxation time is responsible for the quasiparticle multiplication coefficient which in turn 
determines the responsivity and detectivity of the detector and its noise characteristics M. 

The energy relaxation time also serves as a pair-breaking parameter in the supercon- 
ducting density of states which is measured in the tunneling experiment |§. Recently the 
energy relaxation time was measured in cuprate superconductors by studying an electronic 
instability at high vortex velocities in the mixed state 0. 

The electron-electron relaxation time of clean superconductors was calculated for the 
three-dimensional case in Ref. ||] and for the two-dimensional case in Ref. ||. In the last 
work all channels of the electron-electron interaction not only the Coulomb interaction was 
taken into account by using the matrix classification of the interaction channels developed 
earlier in Refs. |7[ and ||. The importance of considering all channels of interaction was 
realized long ago for the problem of gauge invariance in superconductors ||. It was also 
emphasized in Ref. || that interference between different channels of interaction cancels 
divergences in the interaction correction to the superconducting order parameter. 

It is known that in normal impure and low dimensional metals the diffusive motion of 
electrons leads to enhancement of the electron-electron relaxation flO ], [11]. The purpose 
of the present paper is to calculate the electron-electron relaxation time in impure two- 
dimensional superconductors using the formalism of Refs. and ||. Earlier attempt to 
study electron-electron relaxation in impure two-dimensional superconductors |12| took into 
account only the Coulomb electron-electron interaction and therefore ignored the other rel- 
evant interaction channels. As we will show in the present paper including all channels of 
the interaction is very important for the electron-electron relaxation, and leads to results 
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qualitatively different from that of Ref. |12 |. Another difference is that we found impor- 
tant contributions to the scattering relaxation time from gapless collective excitations which 



where mised in Ref. ||12|| . Note also that we derived the quasiparticle energy relaxation time 
from the quantum kinetic equation, not as an imaginary part of the self-energy as in Ref. 

Then we study the superconductor-normal metal two-layer system. Such a system was 
already studied in Ref. || for a clean system. We consider a disordered case in the present 
work and show that recombination relaxation rate is strongly enhanced due to inter-layer 
electron-electron interaction. 



II. ELECTRON-ELECTRON INTERACTION IN MATRIX FORMALISM 

We use the Keldysh diagram technique for nonequilibrium processes in which the electron 
Green's functions, along with the electron-electron interaction potential, the electron self 
energy and the polarization operator are represented by supermatrices 
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The matrix electron Green's function in an impure superconductor in the Nambu represen- 
tation has the form 
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where fj are the Pauli matrices, m is the electron mass, and 



6 « = c fl + _l_) i A^ = Afl+ ' 



2r^ e / V 2r^ e 

(E R ) = (e R ) 2 + (A R ) 2 , (3) 



where A is the energy gap, r is the electron-impurity relaxation time, and 

£ e = (e 2 -A 2 ) 1/2 sgn(4 |e| > A. (4) 

In a spatially uniform system the kinetic components G c and tl c are satisfied the equations 

G C (P) = S(e)[G A (P) - G R (P)}, 
£(P)=S(e)[£ A (P)-£ R (P)}, (5) 

where S(e) = — tanh(e/2T) = 2n(e) — 1, and n(e) is the Fermi distribution function. Similar 
relations hold for matrix interaction potentials and polarization operators: 

y C (Q) = (2N(u) + 1)[V R (Q) - V A (Q)], 
fl c (Q) = (2N(u) + 1)[U R (Q) - n A (Q)], (6) 

where Q = (q, a;) and N(u) is the Bose distribution function. 

For averaging over impurity position it is convenient to introduce the following expres- 
sions 

vf A = — < r 3 G A (P)f t G A (P + Q)f s >, 

1XVT 

vt R = — < T 3 G A (P)T, t G R (P + Q)% >, 

liVT 

where v is the two-spin electron density of states. Calculations give e.g. 

V AA =(±[(1-A)t 3 + Bit 2 ], (8) 



where 



e(e + cu)-A 2 e(e + cu) + A 2 



(e + to)A-eA to A 



B 



(e + uj)A + eA (2e + to) A 
B+ = — = —r- > ( 9 ) 
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{=l + i(Z e+u -Z e )T-Dq 2 T, 
C + = l-i{U. + Qr-Dq 2 r. 



(10) 



Values of rji for different matrices fj are presented in Tab. 1. Note also that the following 
identities hold: A 2 - B 2 = 1 and A\ - B\ = -1. 

Treating the electron-electron interaction in superconductors we use the matrix formalism 
developed in Refs. 7 and 8. The bare vertices for the electron-electron interaction are 
classified in terms of the Pauli matrices. Physical meaning of the corresponding operators 
Oi = ^t$> is the following. Matrix f\ corresponds to the order parameter amplitude A, 
matrix f 2 corresponds to the order parameter phase 0, matrix T3 corresponds to the electron 
density, and the vector matrix kf corresponds to the electric current, the later will not be 
considered in the present paper (see Ref. 5). Note also that each impurity vertex carries the 
matrix f 3 . 

Therefore each interaction vertex operates in both Keldysh and Nambu spaces and has 
the form 7^ n (fi) = 7(7^)^^, where % indicates the component in the Nambu space and 
tensor K^ n stands for the Keldysh space, k is a boson index and m and n are the electron 
indices. In the representation corresponding to Eq. 1 the nonzero components of tensor 



We will omit coefficient l/y/2 in intermediate equations for the impurity renormalized ver- 
tices and restore it in final equations for the polarization operators and the electron self 
energies. 

Impurity averaging leads to the ladder equation for the scalar vertex r(fj) shown in Fig. 
1. Such an equation should be written for each bare matrix fj in the Nambu space. The 
solution of these equations for the vertex r(fj) in the Keldysh-Nambu space is obtained 
following Ref. 13. We start with the equation for the vertex f^^), 




-^22 — Ku — K12 ~ K 2 \ — /-• 



(11) 



< hG A {P)t\ 2 {h)G A {P + Q)f 3 > -f L(r 3 ) = -r 3 . 



(12) 
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The solution of this equation is 



(13) 



Note that renormalized vertex f 22 (f 3 ) has components proportional not only to to matrix 
f 3 but also to matrix f 2 . The other vertices f 22 (fj) and f 21 (fj) are presented in Tab. 2. The 
vertices with the other Keldysh indices are obtained from the equations 



12\'i. 



= (f^fOr, f \ 2 (n) = S(e)[fi 2 (rO - f? 2 (fi)], 
^i(n) = -S(e + io)[^ 2 (n)-Tl 1 (r l )], 
f? 1 (f i )=5(e)f1 1 (f i )-5(c + o;)f? 2 (f i ) 
-[S(e)~ S(e + u)](f 1 22 (T i )y. 



(14) 



Note also that index structure of the renormalized vertices in the Keldysh space is dif- 
ferent from the index structure of the bare vertex described by Eq. (11). It is important 
that the complex conjugate in Eq. (14) operates only in the Keldysh space, thus any i in 
Tab. 2 originating from fj matrix algebra are not affected by the complex conjugate, e. g. 
(compare with Eq. (13)) 



(r 22 (r 3 )y 



(15) 



2 1-C* 

The effective screened electron-electron interaction in a superconductor according to Refs. 
6-8 is shown in Fig. 1. The solution of this equation in 3 x 3 matrix form is 
/-(2/A + Iln)- 1 \ 

[(\/o)- 1 -n 33 ]/© n 23 /p 
n 32 /i? -(2/A + n 22 )/T?y 



V = 



\ 



(16) 



where 



v = -(2/a + n 22 )[(\/ )- 1 - n 33 ] - n 23 n 32 , 



(17) 



and A is the BCS coupling constant (A > 0), Vq — 2ne 2 /q is the nonscreened two-dimensional 
Coulomb potential. 
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The polarization operators renormalized by impurities are expressed through the vertex 
T by the equation 

de 



nS(Q) = -\™t J ^Sp(f S 7--f 3 r? 1 (fi)). 



(18) 



Using Eq. (14) for r^ 1 (f i ) and Tab. 2 we find the polarization operators, 

-S{e + u) S(e) 



K( Q ) = ^ / &((i - .u- su / ) / ( ' ) - (1 + ^ 



n^(g) = ^/&((i-A) 



5(e + cu) -S(e) 



i-C* 



n^(Q) = -!/-— de + 



S(e + w) - 5(e) 

i-C* 



i - C+ i - C 

5(e + cu) S{e) 

i-c; ~i-c+. 

S(e + u;) 5(e) 



I -a 1 - ^ 



5(e) 



(19) 



(20) 



(21) 



(22) 



C* 1-C+ i-c + - 

To calculate the electron relaxation we need the imaginary part of the potentials (the po- 
larization operators) in the quasiparticle representation. Making a transformation from the 
electronic representation to the quasiparticle representation we use Eq. (4). As a result we 
separate out the processes of scattering and recombination of quasiparticles in Eqs. (18)- 
(21). For the imaginary part of the polarization operators we have 



lmH-i{Q)scatt = ^ J A de[S(e + uj) - S(e)]Cu(q,e,u}), 
lining)™ = V -Q(u - 2A) jT~ A de[S(e - u) - S(e)}C u (q, e, -u), 



(23) 
(24) 



where Cu are 



C n (q,e,u) = (l - 



e(e + uj) + A 2 


\ Dq 2 


£e| \£e+u>\ 


)(\M-\Ce\) 2 + W) 2 


e(e + co) + A 2 \ 


Dq 2 


MM ) 


^ £+U} \ + \^\y + ( D q 2 y 


e(e + lu) - A 2 


\ Dq 2 


\Ce\ \Ce+uj\ 


)(\M-\Ce\) 2 + (Dq 2 ) 2 



(25) 



- 1 + 



e(e + a;) 



Dq 2 



(26) 



C 33 (q,e,uj) = -(1 + 
+ 1 ! 



Dg 2 



e(e + w) - A 2 

e(e + co)-A 2 \ Dq 2 



MM J(\M + \&) 2 + wr 
For the off-diagonal polarization operator we have 



(27) 



v 



n 32 (<5)scatt — — 77^ 
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de 



+- 
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S(e) 
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i\£ e+u \-i\S e \+Dq* 
S(e + to) 



(28) 



n&(Q), 



-6(u;-2A)-Au; 
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5( e ) - 5(6 - u) 



+ 



5(c) 



i\^\-i\Q + Dq 2 
S(e-cu) 



(29) 



i|&_ w | + + Dq 2 -i\£ e - u \ - + Dq 2 . 

As we will see in the next chapter, calculating the electron relaxation time we need the 
imaginary part of the propagators Im(V^(Q)). The imaginary part of the propagators may 
originate from the poles of the propagators which correspond to collective excitations or 
from the imaginary part of the polarization operators lm(H^(Q)), which correspond to real 
processes of scattering and recombination of quasiparticles. We will restrict our calculations 
to low temperatures T << A, where large frequencies u — 2 A << T, are important for both 
recombination and scattering processes and small frequencies uj « T « A are important 
only for the scattering processes. The imaginary part of diagonal polarization operators 
for in these regions are presented in Tab. 3. Real parts of the polarization operators are 
analyzed in Appendix A. 

Now we study in detail each of the matrix elements of V. Following Eq. (A17) the 
imaginary part of the potential Vii for scattering processes and for small arguments uj « A 
and Dq 2 « A is 

1 



lmV£(Q) = -Im 



2/A + Il£(Q) 



I / 1 \ 2 



K Im ,/2- i i m n A(Q ) 53 UJ Imn »«»- < 30 > 

This approximation is justified because there is no singularity in the order parameter ampli- 
tude propagator V n at small frequency and momentum for any finite A, which means that 
fluctuations of the amplitude of the order-parameter are massive. For the recombination 
processes, large frequencies u > 2A are important, and according to Tab. 3 Imllii may be 
neglected. 

As was shown in Appendix A, the propagators V33 and V22 for uj « A and Dq 2 « A 
have the form 

A _ K 7T A/y - ^ 2 

V *W- vq*ADq K -(u-iQ)*' 1 j 

V A (0) = 4A 2 ixADq 2 - u 2 

22W) vq7rADqK-(uj-iO) 2 nADq 2 + uj 2 ' 1 ' 

Thus the imaginary part of the propagators V33 and V22 for small arguments comes from 
the pole corresponding to the phase mode, not from the imaginary part of the polarization 
operators. 

For large frequencies cu ~ 2A we use the following approximation, 

lmV£(Q) = (V -'(,) - Ui(Q) - J^M^L. y 1 K v s ( q )Hmn*(Q), (33) 

Im^(Q) = Im(-§ - K(Q) + „ ^^(Q). (34) 

where Vs(q) = K>/v(q + k) is the statically screened Coulomb potential in the normal state, 
which was presented above for the two-dimensional = 2nve 2 is the screening mo- 

mentum. Such an approximation is justified due to absence of collective excitation in this 
frequency region. 
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III. ELECTRON-ELECTRON RELAXATION 



The kinetic equation for nonequilibrium distribution function in a spatially uniform sys- 
tem is 

^ = / ^™[G A {P)][t c {P) - S{e){t A {P) - ± R {P))\. (35) 

The electron energy relaxation time r e _ e is determined from the equation 

1 d dn(e) 

r e -e(T, e) = ~dn(e) dt 

The electron self-energy is shown in Fig. 3. Using the results of Section 2, we have 

1 2& r dQ 



(36) 



x^{Im\//(g)Re[Pr fi (f? 1 (f i ))] - 2Im^Im[Pr f3 (f? 1 (f 2 ))]}, (37) 

where Pr^ means the component proportional to matrix fj (projection on fj). In Eq. (37) 
summation on repeated indices is implied. Using Tab. 2 and relation 5Tl 1 /5S(e) = — 
(r^)* we present Eq. (37) in the form 

;-^ = ^W(§^ H+ " ( " +£)l 



X 



in<(Q)((| + 



eUu, J (Uu, + 6) 2 + (Dq 2 ) 2 
& e(e + co) + A 2 \ Dq 2 



s J(U.-Q 2 + (Dq 2 ) 2 

, e(e + c)-A 2 \ Dq 2 



+ im^(g)((^ + !ii±^). 

we et e+w / 



eU. J (Uu, + 6) 2 + (Dq 2 ) 2 
& e(e + cu)-A 2 \ Dg 2 



e eU* J(U.-Q 2 + (Dq 2 ) 2 

& e(e + cu)-A 2 \ Dq 2 



+. m v- (Q) ((f-i<i±^). 

we ec e+aJ / 



eU. J (Uu, + 6) 2 + (Dq 2 ) 2 
^ + e (e + uj)-A 2 \ Dq 



e eU. J{U.-Q 2 + {Dq 2 ) 2 

A ujA ( Dq 2 Dq 2 

+ 1 2 "^uMu. + &) 2 + (Dq 2 ) 2 + - Q 2 + W) 2 



(38) 



In order to separate out the processes of scattering and recombination of quasiparticles 
we need to make a transformation from the electronic representation to the quasiparticle 
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representation in Eq. (38). Note that the presence of imaginary factor % in the last term 
in Eq. (38) means that for the contribution of the nondiagonal channels of interaction 
requires the states under the gap to be taken into account according to the equation: ( e+LU = 
i[A 2 — (uo + e) 2 ] 1//2 , |e + uo\ < A. Such states should also be included in equations for the 
polarization operator II23, in Eqs. (28) and (29) only the states above the gap were included. 
The analysis similar to that presented in Ref. 6 for the clean case shows that contribution 
from the nondiagonal channels of interaction may be neglected. 

For electrons on the Fermi surface, e = A, 

I 2 f°° f°° Do 2 

= — / duo dqq[N(uo)+n(uo-A)} ' 

r e _ e (T,e = A) scatt J2A Jo 



uo 



Auo + (Dq 



2\2 



+ - 



7T 



X 



lmV 1 i(Q) scatt 
1 



A \ l / 2 

[lmV£(Q) recom + ImV 3 i(Q) recom ] ( — J 

roo roo ]Jq 2 

Jo duj i dqqmuj) + n ^ + A)] 2Au + ^ + W y 

1 ( A \ 1/2 1 



2A + U ) + [Im^(g) scati + lmV£(Q) scatt ] i^— 



00 
2 /■ '■■ 



(39) 



T e _ e (T, . 



X 



- = A) 

/ recom 

scatt 



roo roo 

— duo dqq[N(uo) + n(oo - A)] -5 — 

7T 5 J2A JO CO 2 — 



Dq 2 



to 



uo 

O roo roo 

-4 / duo J 



7T J J2A 



Aa; + (Dq 2 ) 2 
2A\ 1 /2 /Ax 1 / 2 
—J + [Iml^(g) scati + lmV 3 i(Q) scatt ] (-j 

Dq 2 ^ 

% Wm ) + „(, + a)| ^ T7TI ^ 

, 4 ,„. ,/A\ 1/! 



X 



[ImV£(Q) r 



+ ImV£(Q) recam ]^—) 



(40) 



Eqs. (39) and (40) describe processes of scattering "two into two" and "three into one" 
quasiparticles correspondingly. 

Further calculations will be performed for low temperatures T « A. It may be shown 
that the most important contribution to the recombination time originates from terms 
V 22 {Q) recom 

and V 33 (Q) 

recom l n Eq. (40), 

1 T / /Az/\ 2 \ / 2A\ . . 

Calculating the scattering relaxation time from terms V2 2 (Q) r ecom and Vzz(Q) recom in Eq. 
(38) we use the approximation of Eqs. (33) and (34) and we use Tab. 3 for the imaginary 
parts of the polarization operators. As a result we get 
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T (_ (\v\ 2 \ ( A\ 

exp -- . (42) 



r e - e (T,e = A) scatt 2nWu\ ' \ 2j J "V T. 

As for the contribution to the scattering relaxation time from terms ImV^Q) scatt and 
ImV^Q) scatt in Eq. (38) we note that for small energy transfers, u « T « A the 
imaginary part of the propagators V33 and V22 originates from the poles corresponding to 
the phase mode as seen in Eqs. (31) and (32)). Integrating these poles over the momentum 
q we get 

1 21 /7r\V 2 /T\ 3 / 2 T 2 

r e - e {T,e = A) scat - 2A2) UJ £W ( } 

Note that the main contribution comes from the propagator V22 corresponding to the fluc- 
tuation of the phase of the order parameter. 

The term Vu{Q) scatt does not have poles corresponding to the collective mode, thus using 
Eq. (30) we find that the scattering relaxation time has a power law divergence, 

1 A /7TT\V2 , A x 

exp -- . (44) 



- 



r e _ e (T,e = A) scatt 2k*Dv\u> J "V T. 

This divergence is similar to the logarithmic divergence of the phase or energy relaxation 
times in the normal impure two-dimensional case 14 . According to Ref. 14 the cutoff fre- 
quency ojq is defined by the relaxation time r e _ e , which physically means that the kinetic 
equation cannot be applied for energy transfers less than l/r e _ e , thus the self-consistent 
solution of Eq. (44) is 

1 (7rA 2 T) 1 / 3 / 2A\ 

exp -— . (45) 



r e - e (T,e = A) scatt {2^Du)^ "V ST J' 

The low-frequency singularity in the scattering relaxation time mentioned above is for elec- 
trons exactly at the Fermi surface, e = A. We note that for the electrons above the Fermi 
surface, e > A the singularity in the scattering relaxation time associated with the potential 
V11 (Q) scatt is weaker but it does not disappear. More accurately such a divergence must be 
regularized directly in the physically measurable quantity e.g. the tunneling conductance. 
However it is not necessary because the contribution to the scattering relaxation time from 
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the phase collective mode, Eq. (43), is more important because it does not have a small 
exponential factor such as that presented in Eqs. (42) and (45). 

The appearance of the nonexponential scattering relaxation at low temperature is a 
direct consequence of the gapless phase mode in two dimensions. In three dimensions the 
phase mode have a gap and the main contribution to the scattering relaxation comes from 
the potential Vu. 

1 12AT 1 

r e - e (T,e = A) scaa ~ A^Df 1 







(^) 1/4exp ( 


,4) 





(46) 



where z/ 3 = mpp/n 2 is the three-dimensional density of states. Again after regularization of 
singularity in Eq. (46) we have 

Aexp (-— J. (47) 



r e . e (T,e = A) scatt n(nD)^u 3 "V hT ) 
The recombination relaxation time is obtained similar to Eq. (41), 

1 1 A X / 2 T / /Az/ 3 \ 2 \ / 2A\ 

r„(r, e = A)_ = jiw^l 1 + vt) ) exp {—) • (48) 



IV. TWO-LAYER SUPERCONDUCTOR-NORMAL METAL SYSTEM 

We consider a system of two disordered electron layers with different density of states, 
1/1,2, elastic scattering times, ri )2 , mean free paths, £± t 2, and diffusion coefficients Di 2 - The 
layers are coupled by the Coulomb potentials, there is no superconducting coupling between 
the layers. 

First we consider the screened Coulomb potentials in the normal state. The nonscreened 
Coulomb potentials within the layer, Vb, and between electrons in different planes, Uq, are 

1 2-7TP 2 

V (q) = 2ne 2 /qe, U (q) = exp(-g&), (49) 

e i <? 

where e and e 1 are the dielectric constants of the electron layer and the inter-layer media, b 
is the distance between layers. We assume that e e\ and we absorbed e into e 2 
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(51) 



In all further calculations small momentum transfers are important, thus we assume 
qb « 1, 

Vq-U = Vb(l - exp(-gft)) = 27re 2 6, \/ 2 - [/ 2 = V Aire 2 b. (50) 

In this chapter the lower indices of the potentials and the polarization operators refer to 
the layer, e.g. Vu means the Coulomb potential between electrons in the layer 1, etc. The 
screened potentials are satisfied the equations 

V n U 12 \ ( V U \ fV U \ /U, W V u U 12 
U 21 V 22 ) \U Vo) \Uo V o )\0 U 2 ) \U 21 V 22 
We will use the definitions: Vu = Vi, V 22 = V 2 , and U± 2 — U 2 \ = U . The solution of Eq. (3) 
is 

_U V - (V 2 - U 2 )U 2 _ Vp - (V 2 - U 2 )^ 

U - -p, V l~ p > V 2 ~ p , 

P = (1 - VqU^I - V U 2 ) - U 2 U ± U 2 « 1 - V^Il! + n 2 - 47re 2 MI 1 n 2 ). (52) 

The polarization operators in each layers for q£i << 1 and uTi « 1 are chosen in the form 
corresponding to a normal state, because for recombination processes large frequencies uj > 
2 A are important, while collective excitations in a superconductor exist only for u « A. 

n - <« = (53) 

The potentials are 

TfA (m _ 1 (iu + Diq 2 )^ + D 2 q 2 ) 

U W) v x {D x + D 2 v 2 lv x )q* icu + Dq* ' 1 ' 

vm=u*( Q r +{1+2 ? d } D2q \ (55) 



where 



iu + D 2 q 2 



v 2 ,\ D1D2 



D= 1 + - + 2M 1 , k 2 = 27re 2 z/ 2 . (56) 

^1 / D 1 + D 2 h> 2 /i>i 
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We assume that layer 1 is in the superconducting state and layer 2 is in the normal 
state. We will calculate the recombination relaxation time in the superconducting layer due 
to inter-layer electron-electron interaction U. From Eq. (40) we have 

777, 7T — - ? / dcu dqq\N(uj) + n(u — A)l — - — — — 

r e _ e (T,e = A) recom 7t 3 7 2 a Jo qq[ y J 1 n u 2 - Au + {D iq 2 ) 2 

/An 1 / 2 
Im^ W )(-) . (57) 

For the imaginary part of the iner-layer interaction we use the approximation similar to 
Eq. (33), lmU A (q,u) = \U(q, 0)\ 2 lmU 2 (Q), where U(q,0) corresponds to the static limit of 
Eq. (53). The imaginary part of the polarization operator U 2 corresponds to the electron 
scattering in the normal layer and thus does not have a small exponential factor typical for 
a superconductor. As a result 

r e _ e (T,e = A) recom An 2 ^ {D 1 + D 2 v 2 /v l ) 2 \D D + y/2Dj V \ T J' 1 J 

We see that relaxation rate in a superconductor-normal metal two-layer system is increased 
by an exponential factor in comparison with the recombination rate in a single layer, see 
Eq. (40). 

V. CONCLUSIONS 

We derived the kinetic equation describing the electron relaxation in two-dimensional 
impure superconductors. For electrons at the Fermi surface, e = A, the recombination and 
scattering relaxation time were calculated for low temperatures, T « A. We took into 
account all channels of the electron-electron interaction in the superconductor. 

We found that the recombination relaxation rate comes from the quasiparticle scattering 
(recombination processes) associated with the fluctuations of the electron density and the 
phase of the order parameter, the propagators V33 and V 2 2- The recombination relaxation 
rate has double exponential smallness at low temperatures (see Eq. (41)), associated with 
exponentially small number of available quasiparticles. 
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The scattering relaxation rate has a power law temperature dependence (see Eq. (43)) 
due to singularity in the propagators V22 and V33 associated with the gapless collective mode, 
the phase mode. The contribution to the scattering relaxation rate from the fluctuations of 
the amplitude of the order parameter, Vn has an infrared divergence similar to the phase 



relaxation time in the two-dimensional normal metal [I4"|, however after regularization the 
corresponding contribution to the scattering relaxation rate has a small exponential factor 
and therefore is less important that the contribution from the collective excitations. 

We also shown that in the superconductor-normal metal two-layer system the recombina- 
tion relaxation rate in the superconducting layer due to the inter-layer Coulomb interaction is 
strongly increased at low temperatures T << A by an exponential factor exp(A/T) >> 1 in 
comparison with a single superconducting layer, this fact may be important for constructing 
superconducting radiation detectors . 

The author is grateful to I. L. Aleiner for valuable discussions and A. V. Sergeev for his 
help at the early stage of the work. 



APPENDIX A: 



In this Appendix we obtain equations for the polarization operators for some limiting 
cases and prove some identities for them. Though some of the results were already presented 
in Ref. 8, for Matsubara frequencies and in Ref. 15 for T = the analysis for continious fre- 
quencies has some advantages and helps us to estimate the imaginary parts of the interaction 
propagators. 

We start with the polarization operators for q = 0. 1123(0, u) may be taken directly from 
Eq. (21), 

^ , \ ivuiA r . ... 

n 23 (o,^) = — — j(u), (Ai) 

where 
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J(u) = 



de fS(e + cu)-S(e) S(e + u) + S(eY 



Then we take n 22 (u;,0) from Eq. (19), 



11,(0,.) = - ^^)-Ay ( e + .) + 



1 | e(e + a;)-A 2 \ ,S(e + ccO + ,S(e) 



and transform it using the identities 



1 + e(e + c)-A 2 = fe+, + £) 2 -^ 



u 



2U 



e+uj 



to the form 



Now recalling the BCS self consist ancy equation 

2 , v r^Ste 



we see that 



2 v r i ^i e J 



| + n 22 (o,w) = ^j(w). 



To transform n 33 (0, cu), 



„ . , !//•,[/ e(e + cj) - A 2 \ 5(6 + w) - 5(e) 
1133(0,0;) + !/ = -- / de 1 + ^-r^ 1 7 , — + 



e(e + u/) - A 2 \S(e + u/) + S(e) 



u 



ie+w + £e 



we use the identities 



1 + 



e(e + u;)-A 2 £ e+U! - & ( e + u e\ 2A 



u 



e+ui 



+ 



i _ e(e + u;) - A 2 = £ £+u) + £ £ / € + u _ e\ + 



2A 2 
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and get 



n 33 (0, uj) + u = -^-J dt{ e -±^S{t + u>) - ^S(e) 



e 



t+U) 



Now we recall another identity, 



and get 



(k [i±^S{t + u:)-^S{t) 



z/A 2 

n 33 (o,w) = -— 



-2w, 



J(w). 



(A12) 



(A13) 



(A14) 



We see from Eqs. (A5), (A8), and (A13) that the following identity holds 

2 



n 33 (o, u) (- + n 22 (o, w)) + (n 23 (o, uo)f = o. 



(A15) 



At low temperatures, T << A and for small arguments uo « A and Dg 2 << A, the 
polarization operators are 



A 



+ n n (Q) = -v 



cu 2 zV n 

1 h 7T— — 

12A 2 8A 2 



(A16) 



A 



+ n 22 (Q) 



a; 



+ 7T 



2 n 



4A 2 4A 2 



(A17) 
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n 23 (Q) = , n 33 (Q) = -v 

Thus for small arguments the following relation holds 



1 + 



to 



2 n 



6A 2 



n 33 (g) ^- + n 22 (g) J + (n 23 (Q)) 2 ~ g 2 . 

Using Eq. (A16) we have for the propagator Vn 



V 11 (Q) = -(2/\ + U 11 (Q))- 1 = u 1 



CO 7T Do 

+ 



2\ -1 



(A18) 



(A19) 



(A20) 



12A 2 8 A 2 

We see that Vn is not singular which means that fluctuations of the amplitude of the 
order parameter are massive, thus the imaginary part of the propagator Vu originates from 
Imn n (Q). 



The screened Coulomb potential is presented in the form 



V £(Q) = W (A21) 

33 w; l-V (g)lI*(Q)' 1 ; 



where 



ftA (« - «> + 20ik = - %U>Ft- «»)■ • <A22) 

The propagator V22 may be written in a similar way 

V&Q) = -, ' a^i. (A23) 

The poles in the propagators V33 and V22 correspond to the collective excitation, the phase 
mode, which in two dimensions is Jl5] given by equation uj 2 = nADnq. 

In quasi-one-dimensional superconductors the nonscreened Coulomb potential is V (q) = 
2e 2 ln(l/ga), qa << 1, where a is a cross-sectional size, and the density of states is v\ = 



l/ftVp, thus the spectrum of the pase mode is |L8 

2e 2 



uj 2 = nADq 2 ln(l/ga) - 1 (A24) 

HVp 



APPENDIX B: 

It is interesting to see how the spectrum of the phase mode changes in different two-layer 
systems. First we consider a system of two identical impure superconducting planes coupled 
by the Coulomb interaction, but no Josephson coupling between the planes, thus the order 
parameters are independent in each planes. For the screened Coulomb potential in each 
layer we have from Eq. (52) 

v(0) = Vo(q)-[Vo 2 (q)-U^q)MQ) , . 

[l-V Q {q)n(Q)]*-[U {q)IL{Q)] 2 ' 
To avoid confusion we dropped lower indices in the Coulomb potential V. The spectrum 
of the phase modes is defined by equations 1 — (V + U )Il = and 1 — (V — U )tl = 0. 
the solution of these equations for qd « 1 is uj + = ftirADnq) 1 / 2 and = (27iADnd) 1 ^ 2 q. 
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These new phase modes are similar to in-phase and out-of-phase plasmons in symmetric 
two-layer clean normal metal system, see |0J and [ |Hjj . 

In a system of two coupled quasi-one- dimensional disordered superconductors the spec- 
trum of the phase modes in the long- wave limit qd « 1 is 

ul = nADq 2 (—\n(l/qa) - lV uo 2 _ = nADq 2 (— \n(d/a) - 1 ). (B2) 

\7TV F ) \TTVf J 

Now we consider a two-layer superconductor-normal metal disordered system. The po- 
larization operators in the superconducting and normal layers according to Eqs. (A22) and 
(53) are 

n«Q) = -v. *f D f (B3) 

nAD s q 2 — (u — i0) 2 



itf (Q) = -"n . 71 ^ 2 2 . (B4) 

The spectrum of collective excitations is determined from the equation 

1 - V (ill + n 2 - 47re 2 MIin 2 ) = 0. (B5) 



For small momenta q << K s , k s = 1^e 2 v s Eq. (B4) leads to 

to 2 = nAD s q 2 A + iuTrA^, A = 1 + — + 2n s b. (B6) 
The solution of Eq. (B5) is a phase mode with small damping 

u = ^AD s A) l/2 q + ^A^ ) 

aKdJ << — <<l < B7) 

The last inequality is satisfied provided A » 1 and D s « D n . This result was indepen- 
dently obtained in p9| . If the opposite inequality is valid, the solution of Eq. (B6) is a 



diffusion mode, 



2 < D s q 2 
iuj + D n q A = 0, — << mm 



i 



(B8) 
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FIGURES 

FIG. 1. r(fj) is the impurity renormalized scalar vertex in the ladder approximation, Vij is the 
screened electron-electron interaction, £ is the electron self-energy in an impure superconductor. 
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TABLES 

TABLE I. Summary of the results for rjf A and rjf R defined by equation Eq. (7). Quantities 
A, A + , B, and B + are defined by Eq. (9). 



AA 



AR 



no 
m 

m 



%[(l-A+)fb-B+fi] 



%[(1 + A + )t 1 + B + t ] 
£±[(l + A)fr + iBt$] 
%[(l-A)f 3 + iBf 2 ] 



§[(l + A+)fb + B+fi] 
§[(l-A + )fi-B+fb] 
§[(l-A)fi2-»Bfa] 
|[(l + A)f 3 -i5f 2 ] 



TABLE II. Summary of the results for impurity renormalized vertices r 22 (fj) and ^(fj). 



TO 



T2 



T3 



r£ 2 (n) 
r|i(n) 
rL(r 2 ) 

rL(r 3 ) 



C+ 

o 







1 + ^T^7T 












Ci 

2T 



4 



1 + C+ l + A 

Z 2i-C 
*^i-C+ 





1 C i + 4 
C+ l-A 



1 + 



TABLE III. Summary of the results for the polarization operators in the regions of the param- 
eters important for the electron relaxation at low temperatures T « A. 



ImlL 



lo - 2A << T « A 



oj « T « A 



{nuT) l / 2 Dq 2 exp(-A/f) 
" u 2 + 2Aw + (Dg 2 ) 2 
../ttT V /2 (2A + u)Dq* exp(-A/T) 
" ' j oj 2 + 2Atu + (L»g 2 ) 2 
-KUjDq 2 
U uj 2 - 2Aw + (Dg 2 ) 2 

« 



(ImH£ 2 = imn^ 3 ) scatt 
(imn^) scati 



a;^exp(-A/r) 
2Au + (Dg 2 ) 2 
vr^V /2 2A J Dg 2 exp(-A/T) 
2Aoj + (Dq 2 ) 2 



(imn^ 2 = i m n A 



33 jrec 



(imn^) r 
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